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Abstract 

We consider a typical problem in Mean Field Games: the congestion case, where in 
the cost that agents optimize there is a penalization for passing through zones with high 
density of agents, in a deterministic framework. This equilibrium problem is known to be 
equivalent to the optimization of a global functional including an L'' norm of the density. 
The question arises as to produce a similar model replacing the penalization with an 
L°° constraint, but the simplest approaches do not give meaningful definitions. Taking 
into account recent works about crowd motion, where the density constraint p < 1 was 
treated in terms of projections of the velocity field onto the set of admissible velocity 
(with a constraint on the divergence) and a pressure field was introduced, we propose 
a definition and write a system of PDEs including the usual Hamilton- Jacobi equation 
coupled with the continuity equation. For this system, we analyze an example and propose 
some open problems. 

1 Introduction 

The goal of this paper, as well as that of the talk that the author gave at the workshop in 
Rome, is to present an idea to define a Mean Field Game (and to write down the corresponding 
coupled system of PDEs) when the moving density is subject to a congestion constraint, i.e. 
it cannot go beyond a fixed value standing for the maximal compression of the agents. We 
will see that this is not completely evident, but recalls what is done in some papers on "hard" 
congestion for crowd motion. In the same papers, the lack of strategical issues in the movement 
of the agents was perceived as a drawback, which motivates the definition of this new model. 

In the next section we will quickly review the existing models for non-strategic crowd motion 
and for deterministic Mean Field Games with density penalizations. Later on, in Section 
3, we will explain how to glue these ideas so as to put constraints instead of penalizations 
and write down an MFG system. These systems usually couple a continuity equation and a 
Hamilton- Jacobi equation in the unknowns p (density) and tp (value function) , but in this case 
a third unknown will appear, the pressure p, and the system will be complemented with some 
conditions on the pressure. Unfortunately, no result is available right now on this system, and 
Section 4 will be devoted to examples and open questions, in particular those arising from the 
comparison with the penalized models. 
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2 About the existing models for collective motion under 
congestion effects 

We review in this section the models that are the closest to the one we look for. First we 
describe a model for crowd motion where the key point is the constraint p < 1, but no 
strategical issue or equilibrium condition is present; then we pass to a sketchy description of 
Mean Field Games with penalization of high densities. We will conclude with a naive, and 
wrong, idea to replace this penalization with a constraint. 

2.1 Crowd motion with density constraints 

Many models for pedestrian motion are based on the idea that individual slow down when the 
density of the crowd is too high around them. This correspond to saying that their speed is 
(or is influenced by) a decreasing function of the density p. This is what is sometimes called 
"soft congestion". 

A more drastic viewpoint is the one presented by Maury and Venel in [12] , where the idea 
is that particles can move as they want as far as they arc not too dense, but, if a density 
constraint is saturated, then their spontaneous velocity field u could be impossible and would 
be abruptly modified into another field v, less concentrating and typically slower. The main 
assumption of such a model is that v is the projection of u onto the cone of admissible velocities, 
i.e. those that infinitesimally preserve the constraint. Maury and Venel were concerned with 
the discrete (microscopic) case, where individuals are represented by small disks. The density 
constraint is interpreted as a non-superposition constraint: particles cannot overlap, but as 
soon as they are not in touch their motion is unconstrained. When they touch, the set of 
admissible velocities restricts to those that increase the distance of every pair of particles in 
contact. 

Later on, a continuous (macroscopic) model has also been established in [10] (see also 
|11| for the latest developments and a micro- macro comparison). 

• The particles population is described by a probability density p £ Vifl); 

• the non-overlapping constraint is replaced by the condition p E K ~ {p <E 'P(ri) : p < 1} 
(where p denotes, by abuse of notation at the same time the probability and its density, 
since anyway we only consider absolutely continuous measures); 

• for every time t, we consider ut : O — >■ M"* a vector field, possibly depending on p; 

• for every density p we have a set of admissible velocities, characterized by the sign of 
the divergence on the saturated region {p = 1}: adm{p) = {u : 17 — > M'' : V • w > 



• we consider the projection P, which is either the projection in L^{C'^) or in L^{p) (this 
will turn out to be the same, since the only relevant zone is {p ~ 1}); 

• wc solve the equation 



on{p=l}}; 




(1) 
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Equation ([IJ is motivated by the fact that the equation satisfied by the evolution of a 
density p when each particle follows the velocity field v is exactly the continuity equation 
dtPt + V • [ptVt) = (with w • n = on 90, so that the density does not exit 51). Here we only 
insert the fact that v is the projection of u. 

The main difficulty is the fact that the vector field Vt ~ Padm{pt)[ut\ is nor regular (since 
it is obtained as an projection, and may only be expected to be a priori), neither it 
depends regularly on p. 

Anyway, we need to make a bit more precise the definitions above. First, we stress that 
the continuity equation dtp + V ■ (pv) = is to be intended in the weak sense: 



for every ip £ C^(r2), 1 1-^ f ipdpt is Lipschitz and, for a.e. t -y ( i^dpt = [ V?/; • vtdpt- 

Jn dt 

(2) 

Moreover, also the divergence of vector fields which are only supposed to be is to be defined 
in a weak sense, and to do that it is more convenient to give a better description of adm(p) 
by duality: 

adm{p) = l^v e L'^ip) : JvVp<0 Vp E H^{n) : p > 0, p{l - p) = . 

In this way we characterize v = Padm(p)[^i'] through 

u = V + Vp, V € adm{p), J v ■ Vp = 0, 

p e press{p) := {p G H^{^), p>0, p{l - p) = 0}, 

where press{p) is the space of functions p used as test functions in the dual definition of 
adm{p). They play the role of pressures affecting the movement. The two cones \7press{p) 
and adm{p) arc orthogonal cones and this allows for an orthogonal decomposition u = v + Vp. 
This also gives the alternative expression of Equation ([T]), i.e. 

dtp + V ■ {p{u~\/p)) ^0; p>0, p{l- p)=0. (3) 

We stress the fact that the orthogonality condition J v ■ \7p = is no more necessary and is 
a consequence of the continuity equation. This may be informally seen in the following way. 
Fix a time to and notice that 

Ptodpta = J Pto> J Ptodpt 

which means that the function t H- J ptgdpt is maximal at i = to- By differentiating and using 
©, we get 

= y Vpto • vtodpu, = J Vpto ■ vto 

(this proof is only formal because nothing guarantees that t i—> J ptgdpt is diffcrcntiable at 
t — to, but this can be fixed and it is possible to obtain J Vpt ■ vt = for a.e. t). 

As wc said, this PDE is particularly difficult to analyze because of the lack of regularity 
of V, and the main tools to prove at least existence of a solution lie in the use of Optimal 
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Transport and Wasserstein distances. We will not give here details on these topics, which are 
quite well known now and are not crucial for the subsequent parts of the paper. The interested 
reader may refer to [5J (THl [H] . We will only describe here how Optimal Transport enters the 
game in the case where the spontaneous velocity u is given by the gradient of a given function 
that agents want to minimize: u = — VZ?. 

Actually, it happens that in this case, Equation ([T]), may be seen as the gradient flow of 
a functional F with respect to the Wasserstein distance W2 (which is the distance on the set 
7^(ri) induced by the minimal value of the optimal transport problem for the quadratic cost 
\x — yp). The functional that we need to consider is 



F{p) = 



f D{x)dp p e K, 
+00 li p <^ K, 



where, again, K = {p G 'P(f^) : p < 1}. To explain what a gradient flow is, let us only say 
that it is the equation that rules an evolution where po is given, and pt moves by selecting at 
every time the direction that lets the quantity F[pt) decrease as fast as possible. In particular 
this may be seen if one discretizes in time the problem, getting an iterative method. We fix a 
time step r > 0, take p^ ~ po and at every step we find a new density through 



pU, e argminF(p) + ^ argmin / Dix)dp ' 



2t 



Then, by doing it for smaller and smaller time steps r and letting t — >■ 0, one recovers a curve 
of probabilities (continuous for the W2 distance) and can prove that it is a solution of ([T]). 
This follows the abstract ideas of De Giorgi and Ambrosio, see [6l [1], which are meant for 
general metric spaces (see also [2]). When they are applied to the case of the Wasserstein 
distance they allow to obtain evolutionary PDEs, as it has been done first in for the Heat 
and Fokkcr-Plank equations. 

An interesting point here is that the constraint p G K may be seen as a limit of 
penalizations as m — > +00. Indeed, one has 

hm l/p". lim hp\\-^.-l' if|IHI-<M-pe/^ 
m^+oo m J m 1 +00 otherwise. 

This means that the functional F may be approximated, as to —> 00 with the functionals Fm 
defined by Fm{p) = J D{x)p{x)dx + J p{x)"'dx. The precise meaning of this approximation 
is to be read in terms of F— convergence (see [5]), and the applications of F— convergence to 
the corresponding gradient flows are presented in [15]. This penalized counterpart will be 
the starting point for the next subsection. Here we only stress that the Gradient flow of the 
functional Fm gives the well known porous media equations (see |14j ) 

77? — 1 

dtp-\7- {p\/D) _A(p™) = 0. 

2.2 Mean Field Games with congestion and density penalization 

Let us give a quick and informal presentation of the equations of Mean Field Games under 
congestion effects. As we said in the introduction, we will only consider the deterministic case 
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(i.e. there will be no stochastic part in the evolution of the agents and we will use standard 
optimal control theory and not stochastic control). The reader may refer to IHUH] and to the 
rest of this special issue; the presentation of the subject that is provided here does not want 
to be exhaustive but only to address the features of the problem that we need in the sequel. 

Suppose that a population of agents may evolve in time, each agent following trajectories 
of the controlled equation 

y'{t) = f{t,y{y),a{t)), t^[0,T] (4) 

a : [0,T] R" being a control that every agent may choose. At each time the goal of each 
agent is to maximize the payoff 

-[ + ds + <i>(y(r)), (5) 

where g is a given increasing function. This means that a is the effort that every agent makes 
to move in the desired direction, and he pays for it (for simplicity, we take a quadratic cost), 
that its position depends on a through the equation Q , and that he tries to optimize the final 
payoff $ but he also pays for the densities of the regions he passes by. In particular agents 
would like to avoid crowded areas. At this first step the density pit, x) is supposed to be a 
given function. Yet, the MFG problem is an equilibrium problem: given the initial density pQ, 
find a time-dependent family of densities pt such that, when every agent selects its optimal 
trajectory so as to optimize ([5]), the density realized by these optimal trajectories at time t 
are exactly pt- 

One can study the optimal control problem given by (|4]) and ([SJ by means of its value 
function 

^{Ux)^^nY^-j^ i^^^+g{p{s,y{s)))^d.s + ^y{T)) : y{t) = x, y' {s) = f{s,y{s),a{s)). 

It is well-known from optimal control theory that ip satisfies the Hamilton- Jacobi-Bcllmann 
equation (in the viscosity sense) 

dtifiit, x) + H{t, X, SIip{t, x)) ^ 0, ip{T, x) = $(a;), 

where the Hamiltonian H is defined through 

lap 

H{t, X, = sup e • fit, x,a)-^-^- g{pt{x)). (6) 

Moreover, it is also well-known that in the control problem, for every {t,x), the optimal 
choice of a{t) so as to optimize the criterion ([5]) starting from x at time t is the control a 
which maximizes in the definition of H given in ([6|) for ^ = yLp{t,x), i.e. which maximizes 

Vv{t,x)-f{t,x,a)-\-^. 

This gives a system of two coupled equations, since Lp solves an HJB equation where p 
appears in the Hamiltonian, and p evolves according to a continuity equation dtp-\-V ■ (pv) = 0, 
where the vector field v{t,x) is given by f(t,x,at.x) and the dependence of at,x on {t,x) 
passes through \/ip(t,x). To be more precise, we can give an explicit example, in the case 
f(t,x,a) = a (this case corresponds to a classical calculus of variations problem, but we 
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will still use the optimal control formalism since we will need to switch back to true control 
problems in a while). In such a case we can compute explicitly the Hamiltonian H{t, x,^) thus 

getting H{t,x,^) ~ ^ g{pt{x)) and the optimal a in the definition of H{t,x,£_) turns out 

to be exactly ^. Hence, the vector field to be put in the continuity equation is exactly Vip. 
This gives the system 

Idtp + V- (pV^) = 0, (7) 
[ip{T,x) = <^{x), p{0,x) = poix), 

where the first equation is satisfied in the viscosity sense and the second in the distributional 
sense. 

In this case, it is also known that a solution of this system (i.e. an equilibrium in the Mean 
Field Game) may be obtained by minimizing a suitable global functional (obviously, up to a 
change of sign, we could also express it as a maximization). Actually, one can solve 

min / / (^\a{t,x)\'^p{t,x) + G{p{t,x))) dxdt - [ ^{x)p{T,x)dx (8) 
Jo Jn\^ J Jn 

among solutions (p, a) of the continuity equation V- (pa) = with initial datum p(0, x) = 
Po{x). When the function G is chosen as the primitive G oi g (i.e. G' ~ g, and in particular G 
is convex) , it happens that the minimizers of this global functional are equilibria in the sense 
explained above (the trajectories that the agents chose according to these densities are such 
that they exactly realize these densities). These functionals recall the functional proposed 
by Benamou and Brenier [3] to give a dynamical formulation of optimal transport, or those 
studied by Buttazzo, Jimenez and Oudet for density congested dynamics in [4]. 

It is clear that we can choose g{p) = p™^^ (for m > 1), which gives G{p) = ^P™', and the 
minimization problem in tends, as m — > +oo to 

min [ [^\a{t,x)fp{t,x)dxdt~[ <i>{x)p{T, x)dx, (9) 
Jq Jn 2 Jn 

the minimization being, again, performed among solutions (p, a) of the continuity equation 
dtp + V • (pa) = with initial datum p(0, x) = po{x) and satisfying the constraint p{t, a;) < 1 
a.e. In particular, the energy is not finite unless po satisfies the constraint po < 1 as well. By 
the way, this problem is equivalent, at least in a convex domain fl, to finding pj^ < 1 so as to 
minimize TW^ (po, Pt) — J ^dpx and then choosing pt as a constant speed geodesic connecting 
Po to Pt- This comes from the fact that the optimal curve connecting po to px and minimizing 
J J ^\a\^ pdxdt is the geodesic between these two measures in the Wasserstein space, and the 
set K of measures with density bounded by 1 is known to be gcodesically convex (i.e., if 
two measures belong to K, the geodesic between them stays in K as well). This concept of 
convexity and the proof that K is gcodesically convex come from [13] . 

The question arises naturally of how to define a Mean Field Game equilibrium problem 
under the density constraint p < 1 instead of the density penalization +p'"~^. 

It happens that this question is not trivial at all. Let us try the simplest definition. We 
say that a curve of densities pt is an equilibrium if 

• it satisfies the constraint ||pt||L°c < 1 for all times t, 
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• when the agents choose their optimal trajectory, solving 

max -j^ + ■■ y'(^) = "(^)' = 2; 

and respecting the constraint ps{y{s)) < 1, the density they realize at time t is again pt- 

Yet, this definition does not make any sense (and would lead in general to non-existence of 
such an equilibrium) because of this fact: this game is a non-atomic game, where the influence 
of every individual is negligible; thus, every agent facing densities pt < 1 can choose with no 
restrictions at all its trajectory, since he knows that anyway simply adding its own presence 
somewhere will not really change the density and violate the constraint. This means that the 
trajectories that are selected by every agent do not take into account the constraint p < I and 
hence there is no reason for the densities realized by these choices to satisfy p < 1. 



3 The MFG system with density constraint 

The main idea of this section, and of this "modest proposal", is the following: use the pressure! 
Express the effects of the density constraints through the existence of a non-zero pressure field. 
This requires changing the definition of the equilibrium. 

• We describe the situation through a pair (p, a) where pt stands for the density of the 
agents at time t and a{t, x) for the effort that the agent located at x at time t makes to 
control his movement. We require that (p, a) satisfies 

dtPt + ^ - {pt {Padni(pt)[at\)) = 0, 

which means that p is advected by the projection of this effort vector field. 

• Considering the projection of at onto adm{pt) a pressure pt appears such that 

Padm(pt)[at] =at~ Vpt, Pt G press{pt), J Padrn{pt)[at] ■ "^Pt = 

(but, again, the orthogonality condition is included in the continuity equation). 

• Every agent wants to optimize his own control problem where the state equation is 
influenced by the pressure: 

max ^^ds + <i>(y(T)), : y'{s) = a{s) - Vp,(y(s)), y{t) = x. 

• The configuration (p, a) is said to be an equilibrium if the original effort field a coincides 
with the optimal effort in this control problem and if the original densities pt coincide 
with the densities realized at time t according to these optimal trajectories. 

In order to write down the equations that correspond to such an equilibrium we first need 
to compute the Hamiltonian of this control problem: 

lap l£P 

F(i,x,e) = sup e • (a - Vp(t, x)) ~ ^ = ^-^.Vp(i,a;), 
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the optimal a in this maximization being exactly ^. This means that the optimal effort at 
{t, x) will be 'V(p{t, x) and that the vector field appearing in the continuity equation will be 
Vip — Vp. Wc finally get to the MFG system 

{dtifi + ^ - • Vp = 0, 
^{T,x)^^x), p{0,x)^po{x). 

It is important to notice that in this case the equilibrium may only be defined in terms 
of a pair (p, a), and that the density only is not sufficient to describe the configuration. The 
reason lies in the fact that the pressure at time t does not depend on the density pt only, but 
on the vector field that we project on adm{pt). This was not the case in the gradient fiow 
framework of crowd motion in |10| , since the vector field to be projected was itself a function 
of p. This will be clearer in the next example. 



4 An example and some open questions 
4.1 An example where nothing moves 

Even if the example presented in this section is very simple, its resolution will be quite sketchy 
and informal. Consider the following situation: given the function $, suppose that there exists 
£ £ R such that |{$ > £}| — 1. Wc will also suppose $ to be regular enough (in particular, to 
be continuous), so that ^ ~ £ on OA. Let us set A = {$ > £} and take po = Ia^ i-e. P[){x) = 1 
for x G A and pq{x) ~ for x ^ A. 

Let us try to guess the solution of (fTO)) . If instead of a density constraint we had a 
power penalization, the solution would be obtained by minimization of a global functional 
// pa^ + ~// ~/ ^dpT- Let us suppose for a while that a similar principle holds for the 
L°° constraint: we could imagine that the solution is obtained by minimizing JJ pa^ —J ^dpx 
under the constraint p < 1. This problem has an easy solution: not moving at all! Actually, 
Pt = po is a solution since it minimizes the first term (no movement, no kinetic energy) and 
maximizes J ^dp among probability measures with p < 1, since none of them can do better 
than pq, which is already concentrated on the best points for the function $. This suggests 
that the solution of ^TU\\ should be obtained with pt = po- 

Can wc conclude that the pair (p, a) given by pt = po and at — is an equilibrium? the 
answer is for sure not. Why? just because if the effort vector field is identically zero, then its 
projection will also be zero, and the pressure as well. This means that we would have Vpt = 0, 
i.e. the agents could move as if there were no constraints. Hence, every agent would maximize 
his own payoff —J^ _|_ ^(^y(T)), : y(t) = x, y'{s) = a{s), and in general the solution 

is not a ~ 0. 

Yet, even if we did not find an equilibrium so far, we can improve our guess. We can claim 
that the equilibrium configuration is such that everybody wants to move, but nobody manages 
to do it. This means pt = po but a ^ 0. It is indeed possible that the optimal effort a is 
exactly compensated by the pressure effects. Let us look for a solution of ([TO)) of this form. 

We need to impose Vip = Vp on [0, T] x A. This is sufficient to satisfy the continuity 
equation with a static measure pt ~ pO. If we consider that p must vanish where p < 1, this 
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imposes p = on [0,T] x A'^ . Hence we get the (strange) equation satisfied (at least formally) 
by the value function ip: 

dt^+^-^ ^0 on[0,T]x A" and dt^ - = on [0,T] x A. (11) 

Wc need to solve this equation complemented by its final condition ip(T, x) ~ $(a;). Once we 
find if, using Vip ~ Vp on [0, T] x A, wc need to deduce a pressure p and to prove that it is 
actually a pressure, i.e. p > and p{l — p) ~ 0. 

An informal way to solve (|lip is the following. For simplicity, let us suppose ^ = (which 
is always possible up to adding a constant to $). Consider the function $ = — 1$| and define 

^{t,x)^sup-J^ Mfl^ds + |,(y(T)) : y{t)^x, y'{s) = a{s). 

Let us suppose that this value function is regular: it must satisfy 

dt<f + ^-^y- = 0, ^(T, x) = l>(a;), ^ < 0, >p{t, x) = for x e OA. 

The last conditions are easy to check: we see from the definition and from $ < that we have 
0; yet, the constant curve y(s) = x gives (fi{t, x) = 0, provided ^{x) = 0, which is the case 
for X e dA. After that one can define 

= (1 - 2Ia)0, 

which means Lp = (p outside A and p = ~-(p on A. It is possible to check that (at least in a 
formal way) p solves (|lip and p is an admissible pressure. 

Obviously in the analysis of this example wc just provided a solution, with no guarantee 
that it is the unique one. 



4.2 Some open questions 

As we said in the introduction, in this paper we only present System ()10|) , its motivations and 
connections with other problems, but we are not able to prove any result on it. In particular 
existence, uniqueness and well-posedness in general arc open questions. 

The usual conditions guaranteeing uniqueness in MFG are based on sign conditions such 

as 

— / (i^i - v?2)d(pi - P2) > 0, 

to be checked on every pair of solutions (ipi, pi), {ip2, P2)- Yet, they do not seem easy to verify 
here, nor to imply in a simple way uniqueness results, even if some estimates are possible. 

Hence, we ignore for the moment these uniqueness questions, and we only concentrate on 
two questions concerning the connection with the penalization and a possible variational 
principle. 
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First question: can we obtain ([TU]) as a limit of penalizations with m oo? Let 

us be more precise. Consider the System ([7]), choosing g{p) = p"'^^ (which corresponds to a 
penaUzation — J in the global functional that is optimized by the equilibrium) . Obviously, 
in this system there is no pressure, but the question is: is it possible to introduce fictitiously a 
function playing the role of the pressure, thus obtaining a 3-tuple (pm, ^m,Pm) solving System 
([7]), and then wonder if they converge, as to — oo, to a solution [p, ip,p) of (fTO|) ? 

To make an example of what we mean, we can decide to write the value function in ([7]) 
as (f ~ (f — p. Wc will choose later how to define p. If we rewrite the first equation of the 
system we have 

dt^ - dtp + i|V^p - • Vp + i| Vpp - p"'-^ = 0. 

This means that, if we choose p satisfying 

-9tP+i|Vp|2-p"-i=0, 
p(T,x) =0, 

the two equations of the system become exactly the same as those of (jlOp and (p{T, x) = $(.t). 
It is natural to wonder whether the function p that we just defined has some chances to 
converge towards the pressure p of the limit system, and in particular if it will satisfy, at the 
limit, the condition p ~ on {p < 1}. This is not evident, and the only reasonable tool seems 
to interpret the condition —dtp + ^jVp^ — p™^^ = as a Hamilton- Jacobi equation for a 
control problem, thus writing 

P{t,x) = inf (^^M^ +p"-i(y(,s))^ ds : yit) = x, y'{s) = 

Second question: is the equilibrium optimal in some sense? Wc know that for the 
L™ penalization the equilibrium is also obtained if one optimizes the global cost ([8|). One 
could think that letting to — > +oo, the constrained MFG equilibrium should solve 

min j^(^\v{t,x)\^p{t,x)^dxdt- j^^{x)p{T,x)dx (12) 

among solutions (p, v) of the continuity equation 9tp+ V • (pu) = with initial datum p(0, x) = 
Po{x), satisfying the constraint p < 1. First, it is important to notice that here we consider 
the true velocity field v and not the effort field a, and we know that the equilibrium in this 
constrained setting cannot be expressed in terms of {pt,vt) only. This means that, if this 
variational principle is true, then one could use it to identify the velocity field v, but would 
still need to decompose it into Vip — Vp. Yet, it would allow to find the curve of densities pt- 

It is an open question whether this variational principle holds or not. It holds in the 
example of the previous section. To check it in full generality, one possibility would be to use 
the equations of pUj) to check whether the necessary conditions for optimality are satisfied. 

Actually, one can write down by duality some necessary conditions for solving (jl2p under 
the constraint p < 1. In order to apply duality, it is well known, from Benamou-Brenier on 
([3]), that it is convenient to use the variables (p, g) with q = pv instead of (p, w). This leads 
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to the optimization problem 



min 11^-7} / "^dpT + sup - / / dtxdp - / / Vxdq + I Xrdpr - I Xodpo, 



{p,q):0<p<lj J 2p 



X 



where the suprcmum over x stands for the constraint given by the continuity equation. It 
may be dualized by interchanging min and sup and, optimizing first over q, one finds that the 
optimal q is given by q = pVx ^-iid gets to 

sup^min^y" {xt ~ ^)dpT -J Xodpo - J j {dtX + ^I^Xp^ dp. 

This allows (we skip all justifications coming from convex analysis) to write the following 
optimality conditions : if (p, q) is optimal then there exists a function x such that q = pVx 
(and hence v = Vx) and 

• dtX + ^l^xP < on {p = 0} (but this is not restrictive, since x uiay be modified on 
p- negligible sets and still satisfy q = pVx ), 

• 5*x+5lVxl' >0on {p=l}, 

• 9tX+5|VxP=0on{0<p<l}, 

• XT — <& > on {pt ~ 0} (same consideration: it is not restrictive), 

• XT — ^ < on {pT = 1}, 

• XT - = on {0 < Pt < !}■ 



The question of the optimality of the solution of pUj) becomes now a question on the 
conditions satisfied by its velocity field. Since we already know v = V<p — Vp, it is sufficient 
to define x = V ^ P ^^'^ investigate whether the above conditions are satisfied or not. In 
particular we have 

dtX+\\^X? = dtV + - ■ Vp - 9tp+ i|Vpp - -9*^+ i|Vpp, 

and the question is whether this quantity is positive on {p = 1} and zero on {0 < p < 1}. 
This question seems linked to the previous one, where the limit m +00 let appear exactly 
the condition = —dtp + ^jVpp = p™^^ 

The conditions to be verified at time t = T are, instead, easy to check, since (xt — = 
(p{T, x) — p{T, x) — ^{x) = —p{T, x) and we know p > and p = on {0 < p < 1}. 
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